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Role of Equation of State in formation of Black hole 
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We study the physical process of gravitational collapse of a perfect fluid with a linear isentropic 
equation of state p = kp, —1/3 < k < 1. We consider a model with ansatz v'(t,r)/v(t,r) = f'(r) 
that give rise to a family of solutions to Einstein equations with equation of state. The solution so 
obtained lead to homogeneous collapse that evolves from a regular initial data and the positivity 
of energy conditions. The collapse terminates in the formation of black hole. We show that as the 
parameter k —»• 1 , the formation of black hole gets accelerated in time, revealing the significance of 
equation of state in black hole formation. 

PACS numbers: 04.20.Dw, 04.20.Jb, 04.70.-s, 04.70.Bw 


I. INTRODUCTION 

When the star heavier than few solar masses has ex¬ 
hausted its internal nuclear fuel that supplies the out¬ 
ward pressure due to gas and radiation against the inward 
pulling gravitational force, then the star cools down - it 
loses its equilibrium, the unbalanced gravity forces it to 
shrink, and the perpetual gravitational contraction be¬ 
gins. The end state of such contraction develops a space- 
time singularity for a wide range of physically reasonable 
initial data (f]. Concerning such a space-time singular¬ 
ity, Penrose formulated the cosmic censorship conjecture 
(CCC) , it states that ’a singularity of gravitational col¬ 
lapse of a massive star developed from a regular initial 
surface must always be hidden behind the event horizon 
of the gravity ’ 0. This conjecture advocates the forma¬ 
tion of black hole (BH) only as against the formation of 
naked singularity (NS) (wherein the time of formation of 
singularity precedes the epoch of formation of trapped 
surfaces). The CCC is fundamental to the well devel¬ 
oped theory and astrophysical applications of black hole 
physics today. 

Oppenheimer and Snyder studied the spherically sym¬ 
metric model of a homogeneous dust cloud that led to the 
general concept of trapped surfaces and the formation of 
BH 0. The formation of the event horizon takes place 
here, well in advance to the epoch of the formation of the 
spacetime singularity, and hence it is necessarily hidden 
behind the event horizon of the gravity, thus forming the 
BH. This model invigorated the life of BH physics but it 
lacks pressure which is a main constituent in the study 
of gravitational collapse. 

It is known that the physical attributes of the mat¬ 
ter field constituting a star are described by an equation 
of state relating the pressure and density of the matter 
field, and hence it is important to know if the BH would 
form for an assumed equation of state for the collapsing 
cloud. The collapsing massive star evolves to having su¬ 
per dense states of matter close to the end stages of the 
collapse where the physical region has ultra-high densi¬ 
ties, energies and pressures. Such an ultra-high density 
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region of collapse needs to be described by a physically 
realistic equation of state but as of now such an equa¬ 
tion of state is not precisely known, and also whether 
the equation of state would remain unchanged or would 
it actually evolve and change as the collapse develop? , 
these are some of the intriguing questions 0. 

To describe the collapse of a massive star, we can 
choose the equation of state to be linear isentropic or 
polytropic after it loses its equilibrium configuration. 
The gravitational collapse of a perfect fluid with a lin¬ 
ear equation of state is of interest from both theoretical 
as well as numerical relativity perspectives. 

Self-similar perfect fluid collapse models with a linear 
equation of state were considered through numerical sim¬ 
ulations by Ori and Piran 0 and analytically by Joshi 
and Dwivedi Q to show how black holes and naked sin¬ 
gularities develop as collapse final states in this scenario. 
Further, Goswami and Joshi studied the case of an isen¬ 
tropic perfect fluid with a linear equation of state without 
the self-similarity assumption, wherein they showed that 
the occurrence of BH and NS evolving from regular initial 
data depends on the choice of rest of the free functions 
available |7j. 

R. Goswami and P. Joshi have studied a special class 
of perfect fluid collapse models (wherein the mass func¬ 
tion is assumed to be separable in terms of the physical 
radius of the cloud and the time coordinate) which gener¬ 
alizes the homogeneous dust collapse solution in order to 
include non-zero pressures and inhomogeneities into evo¬ 
lution 0. It is shown that a BH is necessarily generated 
as end product of continued gravitational collapse. 

Our motivation also comes from certain other ques¬ 
tions such as, what if the value of k increases in the 
range —1/3 < k < 1 when a BH appears as collapse 
final state for a given spacetime dimension, will the BH 
formed sustain its nature? If, it is so, will the formation 
of BH precede in time as k —> 1 ? We believe answers 
to these and similar issues would be important to un¬ 
derstand better the physical aspects and the role of an 
equation of state in gravitational collapse of a star. 

Therefore, our purpose here is to examine our motiva¬ 
tion for a special class of solutions of the Einstein field 
equations for a spherically symmetric gravitational col¬ 
lapse of perfect fluid to know explicitly how a homoge- 
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neous density profile should behave in the later stages of 
collapse and near the singularity as k —> 1, so that the 
collapse end state would always be a BH. 

In section II, the collapse with linear equation of state 
is studied through the special solution obtained that ter¬ 
minates into BH. The formation of apparent horizon is 
investigated for parameter k —> 1 in section III. The star 
model is completed by studying matching conditions in 
section IV. The conclusions and remarks are specified in 
section V. 


II. HOMOGENEOUS COLLAPSE WITH 
LINEAR EQUATION OF STATE 

In comoving coordinates {t, r, 9, <fi), the general spher¬ 
ically symmetric metric, 

ds 2 = -e 2 ^ r) dt 2 + e 2 ^ r) dr 2 + R 2 (t, r)dfl 2 (1) 

describes spacetime geometry of a collapsing cloud where 
dtd 2 = dd 2 + sin 2 9d(j) 2 is the metric on a two-sphere. 
The stress energy-momentum tensor for the type I matter 
fields for perfect fluids is expressed by 

T” = diag[—p,p,p,p\. 

Herein the quantities p and p are the physical entities rep¬ 
resenting energy density and pressure respectively. The 
weak energy condition is a requirement that for every 
future directed time-like vector field and the stress- 
energy tensor T^ v satisfies the relation T^u^u 1 ' > 0 
which implies p >0; p + p >0. Clouds perfect fluid re¬ 
lation is described through the linear equation of state 


In the continual collapse of the star, we have R < 0, 
it specifies that the physical radius R of the collapsing 
cloud keeps decreasing in time and ultimately, it reaches 
R = 0, and it denotes spacetime singularity, namely the 
shell-focusing singularity at R = 0, where all the matter 
shells collapse to a vanishing physical radius at the epoch 
t = t s . 

The mass function F(t,r ) acts suitably at the regular 
center so that the density remains finite and regular there 
at all times till the occurrence of singular epoch. The 
Misner-Sliarp mass function for the cloud can be written 
in general as, 

F(t,r) = r 3 M(r,v) (8) 

where the function M ( r, v) is regular and continuously 
twice differentiable. On using equation ((HI) in equation 
©, we have 

3 M + r[M. r + M <v t/] M v 
P v 2 (v + rv') k v 2 

We rearrange terms in equation © and express it as 
k rM t r + [(k + 1 )rv' + v)M } v = —3 kM . (10) 

Next to obtain the general solution of equation m 
we consider here the ansatz, 

-=£'(r) (11) 

v 

due to which the equation (HOI) is written in the form 
krM t r + [(k + 1 )rv£’(r) + v]M i v = —3 kM . (12) 

This equation has a general solution of the form, 


p{t,r) = k p(t,r) where k e (—1/3,1]. (2) 


The Einstein field equations for the metric © are written 
as (87 tG = c = 1) Q 


F' 1 F 

P ~ R 2 R’ ~ k R2R 

u, = -YTi M p)] ' 

R'G - 2 RG v' = 0 
F = R [1 — G + H] 


( 3 ) 

( 4 ) 

( 5 ) 

( 6 ) 


where the functions G and H are defined as G(t , r) = 
e~ 2 ^’R' 2 and H(t,r) = e~ 2u R 2 and the arbitrary func¬ 
tion F(t,r) has an interpretation of the mass function 
for the star. On any spacelike hypersurface t = const., 
F(t, r ) determines the total mass of the star in a shell of 
comoving radius r. The weak energy conditions restricts 
F, namely by F(t,r ) > 0, and we have F(t, 0) = 0 to 
preserve the regularity of the model at all the epochs. 

We introduce a new function v(t,r ) by v{t,r) = R/r, 
and using the scalingindependence of the comoving co¬ 
ordinate r, we write |6| 


R(t,r) = r v(t,r). (7) 


e [3(fc+l)«(r)] 

M(r, v) = m 0 -- (13) 

where m 0 is an arbitrary positive constant, and £(r) is 
a continuously differentiable function restrained by no- 
trapped surface condition at the beginning of the collapse 
and by the compatibility condition. The introduction of 
the equation (1111) demands its compatibility with other 
field equations or their subsequent equation (1 25 f) . which 
is discussed in the Appendix. 

M (r, v ) expressed in equation (1131) represents many 
classes of solutions of equation m but only those classes 
are physically realistic which satisfy the energy condi¬ 
tions, which are regular and which give p — > oo as v —> 0. 

Using equation m we write 

e [3(fc+i)£(0)] 

M 0 (v) = M(0, v) = m 0 --, 

M\(v) = -M r (0,u) = 0 (14) 

under the conditions £(0) = const., £'(0) = 0. Here 
M i (r)) = 0 is in accordance with the requirement that the 
energy density has no cusps at the center. The density 
profile for such class of models takes the form, 


. , 3 m 0 e [3(fc+1)J(r) l 

PM = -^(feTT)- 


(15) 
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On integrating equation HID. we obtain 

v(t,r) = e^Q{t) (16) 

where Q{t) is an arbitrary function due to integration. 

We have used the function v(t,r) as a catalyst to find 
solution of field equation ©■ Hereafter, we use (t, r) 
coordinates in understanding the collapse of the dense 
star. The physical quantities in (t, r) take the form 


R(t, 7') = rv{t, r ) = re^Q(t), 

(17) 

F(t,r)=m 0 r 3 e 3 ^Q(t)~ 3k , 

(18) 

pit) = 3 m 0 <5(t) _3(/c+1) . 

(19) 


It is easy to verify that these equations together satisfy 
field equation ©. 

At the dynamical equilibrium event t = U, R(U,r) = 
re z( r ] Q(U) = ri, and 0 < r\ < rb where rf, is the radius 
of the collapsing cloud. The regular density distribution 
at the initial surface takes the form 


In the study of Einsteins field equations with equation 
of state, the system of equations gets closed but still, 
we have introduced equation ED so it needs its com¬ 
patibility with the field equations. It is found that for 
the case b(r) = 0, the function £(r) remains arbitrary 
in satisfying the compatibility condition. Now for such 
physically realistic function £(r), the mass profile of the 
collapsing cloud will be known in r and thereafter, solv¬ 
ing constraint equation ©, we can obtain the function 
Q(t) : thus giving rise to an exact solution of the system 
of equations. 

While for the case b(r) ^ 0, the choice of function £(r) 
is restricted by the condition b(r) = ±b 0 e 2 ^( r ' > where b a 
is a positive constant, thus shrinking the domain of the 
solution set. 

Now, in order to determine Q{t), the field equation © 
can be expressed in the form 


a 2 


F 

R 


G —1 


(25) 


Poiti) = 2>m 0 Q{ti) 3(fe+1) . (20) 

Since po > 0, we have m 0 > 0 with Q(U) = const. 
Further, at some t = tb, r = ri, which is the boundary 
of the cloud where pressure is zero, and where the inte¬ 
rior is matched to a suitable exterior metric, and this is 
discussed in section IV. 

The density profile of the collapsing cloud is homo¬ 
geneous at all the epochs. Now since the particles 
are alike, as the collapse of the star evolves, the den¬ 
sity increases relative to time only on subsequent time- 
hypersurfaces and this feature will be specified by Q(t) 
through Q(t ) < 0. The collapse condition R < 0 becomes 
Q{t) < 0 . 

At the singular epoch t = t s , Q(t) should converge to 
zero and so that density would diverge as t —> t s . So, 
next we aim to find such Q(t) satisfying all the above 
conditions. 

Since the density profile is homogeneous, on integrat¬ 
ing equation (U), we obtain the metric function, 

v(t,r) = a(t). (21) 


where G— 1 = r 2 b(r), and so b(r) basically characterizes 
the energy distribution for the collapsing shells. 

Using equations ED and ED. this can be written as, 


Q(t) = 


= - P AW, 


m 0 e 


2{(r) 


Q(f)U+ 3fe ) 


b(r) 


(26) 


where the negative sign is chosen since for collapse we 
must have Q < 0. Integration of the above equation 
gives us, 


t(r, Q) = ti + 




(27) 


where the variable r is treated as a constant. From equa¬ 
tion ED, we have a regular time function at the center 
of the cloud 


i(0, Q) — ti + 



P «(o) 


dQ 


m 0 e 2 £(°) 

Q(t)( i+3fc) 


+ b(0) 


(28) 


The time for other collapsing shells to arrive at the sin¬ 
gularity can be expressed by 


and using time scaling freedom, we set a(t) = 0. 

We know that the shell-crossing singularity is gravita¬ 
tionally weak and spacetime can be extended through it 
[I]. Therefore, we assume R' > 0, hence from equation 
©. we have 


G(f,r)=d(r) (22) 

where d(r) is another arbitrary continuously differen¬ 
tiable function of r. Further, we write 

d(r) = 1 + r 2 b(r) (23) 


where b(r) is at least twice continuously differentiable. 
The metric takes the form, 

ds 2 = —dt 2 + -- ^ dr 2 -f R 2 (t, r)dfl 2 . (24) 

1 + r 2 b(r) v ' v ’ 


t 8 (r) = t(r , 0) =ti + 



pi(r) 


dQ 


m 0 e 2 €( r ) 

Q( £ )(l + 3fc) 


b(r) 


(29) 


Since the energy density has no cusps at the center, that 
means Mi = 0. The singularity curve then takes the 
form, 


t s (r) =t 0 + rxi(0) + X 2 (0) + £>(r 3 ), (30) 


where Xi(<3) = ~r 
ar 


r=o’ X2(Q) dr 2 


r =0 


and to = i(0, 0) is the time at which the central shell 
becomes singular and which is obtained as 


t{ 0, 0) — ti V 


rQ(ti) 



o 


(31) 
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The time taken by the central shell to reach the singular¬ 
ity should be positive and finite, and hence we have the 
model realistic condition (MRC), namely that, 


m Q e 


2 ?( 0 ) 




+ b(0) 


> 0 i.e. Q(t) (1+3fc) < - 


m 0 e 


2 £( 0 ) 


b(0) 


(32) 


and it must be finite for any k £ (—1/3,1]. 

From consistency of field equations, we have b(0) = 
±6 0 e 2 ^^ 0 - 1 but at the same time Q(t) is a positive real 
valued function, and above inequality makes this possible 
only when b(0) = —b o e 2 ^ 0 \ giving rise to the range of 

Q(t), 0<Q(i)< [m 0 /b 0 ] 1 K 1+sk '> . 

Thus the initial data of mass and density profiles is re¬ 
stricted by the introduction of the equation HID through 
the condition b(r) = —b 0 e 2 ^ r \ 

For b(0) = 0, the above MRC takes the form Q{t ) > 0. 
Now for various choices of £(r), we can study different 
models, for eg. in the case of b(r) = 0 if we choose 
£(r) = 0, metric (l24l) gives us Einstein-deSitter model 
with equation of state while for b 0 = l,£(r) = 0, we have 
a closed Friedman model, and so on. 


Now using b(r) = —b 0 e 2 ^ r \ equation (l27l) takes the 
form 


t{Q) — ti + 



dQ 


Q(i)( 1+3*0 


(33) 


On solving above equation, we have 
2 


t(Q)= U 


. 3(fc+l) ^ . . 3(fc+l) 

Q(U) 2 Hi - Q(t) 2 H 2 


3 v / m0(l + k) 

where Hi= hypergeom([l/2, AT], [AT], b 0 Q(U) K3 /m 0 ), 
H 2 = hypergeom([l/2, AT], [K 2 \,b 0 Q(t) K3 /m 0 ), 
3(fc + 1) {9k+ 5) 


(34) 


AT = 


2(3fc + 1) 


, AT = 


2(3fc + 1) 


and AT = 3k + 1. 


The hypergeometric series mentioned above is convergent 
for \b 0 Q{t) K3 /m 0 \ < 1 and —1/3 < k < 1. The conver¬ 
gence condition on Q{t ) augurs well with the MRC re¬ 
striction 0 < Q{t) < [?7i 0 /& 0 ] 1 ^ 1+3fc ^ ■ We can find Q(t) 
from above equation 


Q{t) = - 

where 


2?n 0 3/2 ATQ(f)- (3fc+1)/2 
[2 m 0 K 2 H 2 + b a H 3 Q{t ) K2 ] 


(35) 


H 3 = hypergeom ([3/2, AT + 1], [A' 2 + 1], b 0 Q K3 /m 0 ) . 


The collapse condition Q{t) < 0 for the dense cloud is 
thus satisfied and as t —> t s , indicates perpetual gravita¬ 
tional collapse of the star. 

The time taken by the central shell to reach the singu¬ 
larity is given by 


to — ti + 


_ , 3(fc + l) 

2Q{t l )^H 1 
3^/m^{l + k) 


(36) 


We observe through the implication of equation (1341) on 
equation (l30l) is that t s {r) = to, indicating that time 


of formation of central singularity (t = t s , r = 0) and 
the non-central singularity (f = f s ,r = r c > 0) in the 
neighbourhood of the center r = 0 is same. Clearly these 
events are simultaneous and it is expected that in such 
scenario the singularity be remain covered behind the 
event horizon. 

In the analysis of gravitational collapse of a dense star, 
we must have the initial configuration to be not trapped. 
Therefore, we must have 

= m 0 r 2 e 2 ^Q(t z )- {1+3k) < 1 , (37) 

R{ti,r) 


this allows for the formation of trapped surfaces during 
the collapse. This constraint indicates, how the choice 
of the initial matter configuration A(T,r) through £(r) 
is related to the initial surface of the collapsing cloud. 
Some restrictions on the choices of function £(r) must be 
made in order, not to have trapped surfaces at the initial 
time. Further, to prevent trapped surfaces at the ini¬ 
tial epoch the velocity of the infalling shells must satisfy 
|i?| > y/\d(r)\. It clearly shows that the initial velocity 
of the infalling shells of the cloud must always be posi¬ 
tive and that the case of equilibrium configuration where 
R = 0 can be taken only at the static boundary of the 
star where pressure is zero Q. 


III. THE APPARENT HORIZON 


When the singularity curve is constant (\i and other 
higher order terms are all vanishing), or would be de¬ 
creasing, then a black hole will necessarily form as the 
collapse final state. For a black hole to come into being 
the trapped surfaces form before the formation of singu¬ 
larity. Thus for a black hole to form we require, [3] 

t a h{r) < to for r > 0 , near r = 0 (38) 

where to is the epoch at which the central shell hits the 
singularity. 

We know since the collapsing shells are simultaneous, 
the end state of collapse is bound to be a black hole or 
to say condition (1551) will be satisfied but the intriguing 
question is that, what is the role of parameter k of equa¬ 
tion of state in the formation of the black hole. Is their 
a certain range of k in which the formation of black hole 
will be accelerated in time? 

The introduction of equation HID has imposed restric¬ 
tion on the function b(r) through compatibility of field 
equations and that now we have b(r) = —b 0 e^ r \ There¬ 
fore, we have only two cases to study namely that space- 
time is bound or marginally bound. The time of oc¬ 
currence of apparent horizon in a bounded spacetime is 
written as 


tbah — tbs 



dQ 


m 0 _ i 

Q(t)( 1+3*0 U 0 


(39) 
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where t bs = t 0 , given by equation fl36j and Q(t a h) = Qah 
obtained from F/R = 1 is specified as 


Qah — 


m 0 r ah 2 e 2f(r “' l) 


(1+3*:) 


(40) 


further we can express 


Rah — 



rah 2 e 2?(r 


3(1 + *:) 
[2(1 + 3*:)) 


(41) 


tbs > tbk for all k £ (-1/3,1]- Further t ba h < hs for 
any r > 0, near the center r = 0 and t b h is a positive 
decreasing time function as k —> 1 . 

Proof: Since R a h remains same for the given mass and 
physical radius of the collapsing star, irrespective of the 
different values of equation of state parameter k, then 
equation (UTl) dictates m 0 and r a h to vary relative to k. 
Therefore using equations (TUI) . C3l) and (PHI) , we have 


On solving equation (l3iH). we have 


dr ah rah ln 

Qti 

Qah 



3(fc+l) 

dk (1 T k)[l + r ah C 

(rah)} 


3(fc+l) atx (,1 -f A.JL 1 "T ' a/ts (I ah Jj 

t b ah= hs - hk where t b k = ir7= 7TTTY ( 42 ) ^ tbk _ ^ b ° R ahr a h e 2 ^ rah ">H 3rah 

3^(1 + k) drah (9k + 5)[l + r ah ?(r ah )]- 1 


(45) 

(46) 


and H 2a h= liypergeom([l/2, K{\, [K 2 ], b 0 Q^/m 0 ). 

It is clear that t b k is a positive quantity for all k £ 
(—1/3,1], and therefore tbah < tbs for any r > 0, near 
the center r = 0. The collapse progresses to culminate 
into the formation of trapped surfaces first and eventu¬ 
ally the singularity forms later, leading to formation of 
BH as a final state of collapse for all k £ (—1/3,1]. 

Now, we study the characteristics of the parameter k 
in the formation of the BH. Let us consider the situation 
where we have known initial mass and physical radius of 
the collapsing star, and this leads to the determination 
of component m 0 of mass or density of the star. And 
since our model is homogeneous in density, therefore this 
density shall be alike at all r at the initial epoch. 

Using equations (PT71) and (PT51) we can write 


to 0 


F(t. u r)Q(U) 3{1+k) 

R 3 {ti,r) 


(43) 


and we find that m 0 decreases as k —> 1 because 0 < 
Q(t) < [ m 0 /b 0 ] 1 ^ 1+3k ) < 1. We can rearrange the ex¬ 
pression for tbk given in equation (l42l) as follows 


; _ 2 R a h H 2ra h 

bk ~ 3(1 + k) 

where H 2rah = hypergeom([l/2, K{\, [K 2 ],z ) 
and *= b Q r ah 2 e 2 «^. 


(44) 


dhk 


2 Rah 


d 

(IT k)-j^j-H 2ra h R 2 rah 


< 0 (47) 


dk 3(1 + k) 2 
where H 3rah = hypergeom([3/2, K l + 1], [I< 2 + l],z) 


and signs are prescribed under the conditions that 0 < 
Qah < Q(ti ) < 1, Q{ti)/Qah > 1 and |z| < 1. 
These physically realistic conditions are possible with 
the appropriate choice of the function £(r) such as 
[1 T r a h£'(r a h)] > 0- Now, we can write 


dtbk dt b k dtbk dr a h 
dk dk dr a h dk 


(48) 


We can obtain dt bs /dk using equations (PHI) and fi3l) . 


dtbs 

dk 


2Q(U) 3(fe+1)/2 

3(1 + k) 2 y/mZ 




< 0 (49) 


From equations (f36| ) and (|42|) , we have tb s > 0 and tbk > 
0, further we can write 


tbs 

tbk 


+ k) ti 

9 tt x-)3(l+fc)/2 
^ J^2ah Q a h 


\ Hl 1 

Q(ti)' 

_ R-2ah _ 

Qah 


3(i+fc) 


> 1 


(50) 

Therefore, from above we conclude that t b k and t bs are 
positive decreasing time functions as k —> 1. 


Now the star with a known radius and mass collapses 
under its own gravitational pull towards the center 
which is solely caused due to its mass, and not because 
of what composition of matter it has. The boundary of 
the horizon of such a star can be known through F = R, 
this gives us a specific radius Rah , independent of the 
equation of state but equation of state can stimulate the 
scenario of formation of trapped surfaces. It is indeed 
possible to testify whether formation of trapped surfaces 
of such a star would accelerate or decelerate in time 
relative to change in parameter of equation of state. In 
view of these aspects the theorem follows: 

Theorem 1. 

Consider t b k = tbk{k,r a h ), r a h depends on k and 
0 < Qah < Q(U) < 1. We prove that both tbk and t b s 
are positive decreasing time functions as k —> 1 and that 


Now tbah = hs - tbk and t bs > tbk > 0, therefore 
tbah < tbs for any r > 0, near the center r = 0. Clearly 
indicating that trapped surfaces are being formed first, 
and the event of the formation of singularity is taking 
place at the later time. Thus black hole forms for all k. 

Further since both t bs and t b k are positive decreasing 
functions as k —> 1 and t bs > t b k- Therefore tbah is a 
positive decreasing function as k —> 1. <C> 


In the marginally bound case that is when b(r) = 0, 
on integrating equation (l26l) . we have 


Q{i) 



2/[3(l+fc)] 


*) 


(51) 


then Q(t) < 0 and as t —> t s , Q(t) — oo. The physical 
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radius and density of the collapsing star are obtained as 


R(t, r ) = re^ 

Pit) = 


+ k)(t s - t) 


2/[3(l+fc)] 


3(1 + k) 2 (t s - t ) 2 


(52) 


The apparent horizon equation using equation (l5lT) is ex¬ 
pressed by 


tah =t s -tk where 

3(l + fc) 

A 2 [m 0 r ah 2 e 2 « r ^)] 
tk — - 


2 Ra 


3(1 + k)^JnT 0 


3(1 + k) 


(53) 


Time taken by the shells to reach the singularity is given 
by equation (l36l) with Hi = 1 ( or by equation (I5T1) ) 


‘• = ‘ (0) = ‘- + 3(TTW« t -> 3,,+W2 - (54> 

Clearly, Theorem 1. holds for marginally bound space- 
time wherein b Q = 0. 

It is evident from Theorem 1. that the equation of state is 
stimulating the formation of apparent horizon of gravity 
to take place at the earlier epoch and further strength¬ 
ening this characteristic as k increases as compared to 
the usual process of formation of trapped surfaces in the 
final stages of collapse of the sufficiently large star, cul¬ 
minating it into the black hole at the earlier time. This 
process is accelerated in time as k —> 1 with the physi¬ 
cally plausible choice of the function £(r). 

To have further insight into the end stages of collapse 
of the star and to analyze the role of parameter fc, we 
consider an example of a Neutron star temporarily in 
equilibrium state, having a mass of 3.24M 0 and physical 
radius of 18.02Ato. This star in due course of time col¬ 
lapses under its gravitational force. The dominant tidal 
force culminates it into the BH with mass of 1.7192Mq 
and the radius shrinks to 5.07A'm where mass F = 2 M 
H3. We have at the horizon, R a h given by equation 
m, thereby r a h is analyzed through numerical solution 
by expanding e 2 ^ r “'*) to the second order, for the phys¬ 
ically plausible choice of the function £(r). The results 
are shown in Table I & II. 

In Table I and II, we have observed that as k —> 1, 
the time of formation of singularity t s decreases in time, 
and that the time of formation of event horizon of gravity 
t a h as well decreases but precedes in time to t s . So the 
trapped surfaces form well in advance in time before the 
event of formation of singularity takes place, culminating 
the final stages of collapse into a black hole. Further the 
formation of BH is accelerated in time as k —> 1 in the 
sense that trapped surfaces are coming into existence at 
the earliest time and thus this property is strengthened 
for increasing parameter of equation of state. 


IV. EXTERIOR SPACE-TIME AND JUNCTION 
CONDITIONS 

The space time in the exterior region denoted as A4 + 
of the collapsing stellar configuration will be filled with 


TABLE I. Role of k is analyzed through the time of formation 
of horizon and that of singularity for initial mass of 3.24 Mq 
and physical radius of 18.02Xm. The collapse end state BH 
has mass of 1.7192 Mq and the Schwarzschild radius is of 
5.07 Km with £(r) = 1 + r 2 /2, U = 0, Q(ti) = 0.37 and 
b 0 = 0.2227 x 10 -3 . 


k 

m, 

o 

T ah 

z 

t s 

tbk tbah — ts tbk 

-0.3 

0.20 

X 

10' 3 

4.216 

0.856 

109.33 

10.65 

98.6848 

-0.2 

0.15 

X 

10 -3 

3.930 

0.587 

58.50 

5.73 

52.7684 

-0.1 

0.11 

X 

10~ 3 

3.719 

0.438 

43.74 

4.47 

39.2577 

0 

0.83 

X 

10~ 4 

3.556 

0.346 

35.79 

3.82 

31.9723 

0.1 

0.61 

X 

10~ 4 

3.427 

0.286 

30.59 

3.36 

27.2268 

0.2 

0.46 

X 

10" 4 

3.322 

0.243 

26.85 

3.03 

23.8244 

0.3 

0.34 

X 

10~ 4 

3.235 

0.213 

23.99 

2.76 

21.2377 

0.4 

0.25 

X 

10~ 4 

3.162 

0.189 

21.73 

2.54 

19.1912 

0.5 

0.19 

X 

10~ 4 

3.099 

0.171 

19.88 

2.35 

17.5245 

0.6 

0.14 

X 

10~ 4 

3.045 

0.157 

18.33 

2.20 

16.1370 

0.7 

0.10 

X 

10” 4 

2.998 

0.145 

17.02 

2.06 

14.9614 

0.8 

0.76 

X 

10~ 5 

2.957 

0.135 

15.89 

1.94 

13.9512 

0.9 

0.56 

X 

10” 5 

2.920 

0.127 

14.90 

1.83 

13.0728 

1.0 

0.42 

X 

10~ 5 

2.887 

0.120 

14.04 

1.74 

12.3014 


TABLE II. Time of formation of singularity and horizon are 
studied for the same data conceived in Table I except that 
b 0 = 0. Herein z = 0. 


k 

m, 

0 

^ ah 

ts 

tk 

tah t s tfe 

-0.3 0.20 

X 

10 -3 

4.216 

23.5566 

4.8286 

18.7280 

-0.2 

0.15 

X 

10~ 3 

3.930 

20.6120 

4.2250 

16.3870 

-0.1 

0.11 

X 

10~ 3 

3.7186 

18.3218 

3.7556 

14.5662 

0 

0.82 

X 

10~ 4 

3.5559 

16.4896 

3.3800 

13.1096 

0.1 

0.61 

X 

10~ 4 

3.4268 

14.9905 

3.0727 

11.9178 

0.2 

0.46 

X 

10~ 4 

3.3219 

13.7413 

2.8167 

10.9247 

0.3 

0.34 

X 

10~ 4 

3.2350 

12.6843 

2.6000 

10.0843 

0.4 

0.25 

X 

10' 4 

3.1619 

11.7783 

2.4143 

9.3640 

0.5 

0.19 

X 

10~ 4 

3.0994 

10.9931 

2.2533 

8.7397 

0.6 

0.14 

X 

10~ 4 

3.0455 

10.3060 

2.1125 

8.1935 

0.7 

0.10 

X 

10~ 4 

2.9984 

9.6998 

1.9882 

7.7115 

0.8 

0.76 

X 

10~ 5 

2.9570 

9.1609 

1.8778 

7.2831 

0.9 

0.57 

X 

10~ 5 

2.9203 

8.6787 

1.7789 

6.8998 

1.0 

0.42 

X 

10~ 5 

2.8876 

8.2448 

1.6900 

6.5548 


radiation flowing outward along radial direction and it is 
appropriately described by the Vaidya metric 

ds\ = - (\ - 2M y V ^ dV 2 - 2 dVdy + yW. (55) 

Let denote the space time in the interior of the 
collapsing star which is separated from the exterior by 
a time like 3 dimensional space time surface £ which 
represents at any instant the boundary separating A4 + 
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from A4 _. The intrinsic metric on E will be 

ds 2 = — dr 2 + 7 Z 2 {r)dfl 2 (56) 

The metric on the interior manifold is described by 


dy 


dV 


v Tt +v h; 



2 M(V) 
y 


E 


Ry/l + ?’ 2 b(r) 


( 66 ) 

The equations m and (1661) determine the mass con¬ 
tained within spherical region in A4+ as 


d/2 

ds 2 _ = -dt 2 + -=-— dr 2 + R 2 (t,r)dtt 2 . (57) 

1 + r z b(r) 



(67) 


The boundary conditions smoothly joining the interior 
and exterior manifold A4- and A4 + across E are stipu¬ 
lated as EMU, 


(ds 2 _) e = (c 7 s+)e = (ds 2 )s 


(58) 


K i: = K ti 


where 


8 2 a 


JC.. = — 77 " 
tj a d?d& 


-n*r 


, „ <9x£ 9 x1 

“ d f 4 d& 


(59) 


(60) 


denote extrinsic curvatures of E in M± respectively 0, 
0 . 

The corresponding normal vectors are 


n„ = 0, 


Rf 


•\/l + r 2 b(r) 




(61) 

The boundary conditions (1581) imply the following rela¬ 
tions 


^ = 1 > »”) = = 2/ 


(62) 


We obtain dV/dr and d 2 V/dr 2 from equations (1621) and 
(1551) , on using them together with equation (KT7D , we find 
that the condition (l65l) leads to the relation 


2 RR 


s 


r 2 b(r) - R 2 



( 68 ) 


Further, on use of the equations ED and m , we can 
write 


(1 + 3 k) m 0 re« (r) Q(t)- (2+3fc) 


(2+3fc) 


F 


E 

R 2 


(69) 


J E 


Finally using equation of state p = kp together with 
equations ED and (LLLLI) • we obtain 


[p QW] e = 0. (70) 

We know at the boundary of the star, the mass and the 
physical radius of the star are fixed numbers, and hence 
at the boundary Q(t )s ^ 0. Hence we must have 


(p)s — 0- 


Thus we could establish a complete star model which ini¬ 
tially match with radiating Vaidya region, and thereof 
with the empty exterior represented by Schwarzscliild 
metric. 


and 




(63) 


V. CONCLUSIONS AND REMARKS 


The extrinsic curvatures Kjj of E are found to have the 
following explicit expressions 


K~ r = 0s 
K+ = 


^ee ~ 


Ke = 


d 2 V 

(dV\ 

dT 2 

Wr J 

R sj 1 + r 2 b 

dy 

dV 




M(V) dV 

y 2 dr 


1 - 


2 M 


J E 


K ii = K ti = 0 for i^j 


(64a) 

(64b) 

(64c) 

(64d) 

(64e) 


In view of Eqs. (I64al) to (I64dl) . the boundary conditions 
ensuring continuity of extrinsic curvatures across E imply 
the following relations 


'd 2 V 


' M(V) 

( dV Y] 

dr 2 

E 

y 2 

\ dr ) _ 


(65) 


Let us summarize the results, firstly we have obtained 
the solution of Type I matter field equations through the 
ansatz introduced in equation ED- Certainly, this has 
led to a special class of solutions with an isentropic equa¬ 
tion of state p = kp that satisfy weak energy conditions 
and evolve as the collapse begins according to the homo¬ 
geneous distribution of matter. 

With the varied choices of function £(r) satisfying 
physically realistic conditions, we have a class of bound 
and marginally bound space-times which can be explored 
further. It is shown that how the choice of initial data of 
mass function and the physical radius through the func¬ 
tion £(r) lead to the formation of BH. 

Studies show that gravitational lensing is an impor¬ 
tant astrophysical tool to observationally test the Cosmic 
Censorship Hypothesis (CCH) 0,0. Relativistic im¬ 
ages of Schwarzscliild black hole lensing is studied by K. 
S. Virbhadra 0], In view of these aspects, we empha¬ 
size that the BH model presented here may serve as an 
example to understand black hole physics in the light of 














































proving or formulating CCH in dynamical gravitational 
collapse. 

The investigation of gravitational collapse with a linear 
equation of state has revealed the role of the parameter k 
in terms of formation of BH and further strengthening it 
by accelerating the formation of trapped surfaces in time, 
in both the bound and marginally bound space-times. 

The parameter value k = 1 depicts the case of stiff 
fluid ( that the equation of state becomes rigid enough ) 
which itself may halt the progress of the collapse at some 
stage 0- Therefore our results are more significant in 


the range of —1/3 < k < 1. 

Though this aspect together with whether the equa¬ 
tion of state would remain unchanged or it would ac¬ 
tually evolve and change as the collapse develops could 
not be established because of complexities. Also the un¬ 
bound case of space-time could not be studied with the 
solution exhibited because of restrictions imposed by the 
compatibility conditions. 
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APPENDIX- COMPATIBILITY CONDITION 


The introduction of the equation CD which is given 
by v'(t,r) = (r), demands its compatibility with 

other field equations or their subsequent equation (1261) . 
So, we consider 

v' = W(t , r, v, v ', v) and v = U(t, r, v, v ', v). (71) 


The condition of compatibility for non-linear partial 
differential equations of order one yields, 


IV, t = U r . _ (72) 

We have W = v(t, r)£'(r) and U = D(t.,r ), above 
equation takes the form 


-y/Dfar) =- D ' {t " r) 

2 VWu) 

»2?(r) 


where D(t,r ) = 


m 0 e 


Q(ty+ 3k 

On simplification, we obtain 

db . . d£ 
- = 2b r -i 
dr dr 


b(r). 


(73) 


(74) 


and solving this equation, we have requisite condition of 
compatibility, 


b(r) = ± b a e 2 ^ r \ 








